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Hidden symmetry in Coulomb interaction is one of the mysterious problems of modern physics. 
Additional conserved quantities associated with extra symmetry govern wide variety of physics 
problems, from planetary motion till fine and hyperfine structures of atomic spectra. In this paper 
we present a simple derivation of hidden symmetry operator in relativistic quantum mechanics for 
the Dirac equation in the Coulomb field. We established that this operator may be reduced to 
the one introduced by Johnson and Lippmann. It is worthwhile to notice that this operator was 
discussed in literature very rarely and so is not known well among physicists and was omitted even 
in the recent textbooks on relativistic quantum mechanics and/or quantum electrodynamics. 



I. INTRODUCTION 

It is well known that the Coulomb problem has 
additional dynamical symmetry in classical and non- 
relativistic quantum mechanics. Classical orbits in the 
Coulomb field are closed due to this symmetry. In 
quantum mechanics this symmetry shows up in the 
form of " accidental" degeneracy of hydrogen atom spec- 
trum. The nature of these phenomena was explained 
by V.Focki, V. BargmanS, W. Paul& and others many 
years ago. The "hidden" symmetry related applications 
have been discussed in several textbooks^^ as well. The 
conserved quantity associated with this symmetry often 
called the Laplace-Runge-Lentz (LRL) vector were dis- 
cussed by many authors (see^, for instance). 

Less information is published about this symmetry in 
the relativistic quantum mechanics, in particular, for 
the Dirac equation. It is known, that well known l- 
degeneracy of hydrogen atom spectrum in non-relativistic 
quantum mechanics, is removed in case of the Dirac equa- 
tion (since LRL vector is no longer conserved in this 
case). However the degeneracy is not lifted completely 
even in this case. The two-fold degeneracy in total an- 
gular momentum still remains and it contradicts with 
existence of well known experimental result - the Lamb 
shift. It is entirely possible that the Lamb shift is also a 
consequence of the same symmetry and it was discussed 
in this content in relatively recent papersSJ^. 

The generalization of LRL vector that involves spin 
degrees of freedom (Dirac's particle) is called Johnson- 
Lippman (JL) operator. The aim of our paper is to de- 
rive the JL operator for relativistic hydrogen atom in the 
Dirac equation as easily as possible. 



II. THE JOHNSON-LIPPMAN OPERATOR (A 
SHORT HISTORY) 

The analogue of LRL vector in the Dirac equation is 
the Johnson-Lippmann (JL) operator which in 1950 was 
published by these authors in the form of brief abstract. 
Because of importance and rarity of this problem, below 



we display the full text of the above mentioned abstract 
as it was presented then. 

"J6. relativistic Kepler Problem, M.H. Johnson 
and B.A. Lippman, Naval Research Laboratory - Besides 
the usual integrals of motion, M and j (in Dirac's nota- 
tion) the relativistic equations for a charge in a Coulomb 
field admit 



A = 3rr 1 - i ( — I (mc 2 ) 1 jpi(H ~ mc 2 p 3 ) (1) 



as another integral of motion. Since A and j anticom- 
mute, the pairs with the same \j\ are degenerate. Thus 
the existence of A establishes the "accidental" degener- 
acy in relativistic Kepler problem just as the existence of 
the axial vector establishes the degeneracy with respect 
to I in the corresponding non-relativistic problem" . 

It is all what is known about this operator in scientific 
press and textbooks. More detailed consideration, in our 
best knowledge, had not been published neither then nor 
after. As far as textbooks concerned, we have found only 
a. brief remark in the form of a short footnote in B.V. 
Berestetskii's booki£. 

Therefore we think that the derivation of JL operator 
and consideration its principal properties is worthy of 
certain interest for students as well as for teachers of 
universities. 

The JL operator were generalized to arbitrary dimen- 
sions in the recent paper— and it was shown that this op- 
erator can be used to construct relativistic supercharge 
in for so-called Witten's superalgebra. As far as com- 
mutativity of JL operator with the Dirac Hamiltonian is 
concerned, it was only mentioned that it can be proved 
by " rather tedious manipulations" . 

We think that from pedagogic point of view it is im- 
portant to show explicitly the simplest way to prove the 
above statement. Below we construct this operator in 
rather simple and very transparent and understandable 
way. 
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III. DIRAC'S A-OPERATOR AND A-ODD 
OPERATORS 

It is well known that the Dirac's operator 

K = /3(S • F+ 1) (2) 

commutes with the Dirac Hamiltonian 

(y. 

H = a ■ p + (im , a = Ze 2 = Za (3) 

r 

where I is the angular momentum vector. Here a and 
/3 are usual Dirac matrices and e 2 = a is is the fine 
structure constant, while E is the electron spin matrix 



£ = pio? = 7°a = 



CT 

a 



(4) 



According Sommerfeld formula the hydrogen atom spec- 
trum is given as follows^: 




(Za) 



(n - \k\ + ^k 2 - (Za) 2 ) 



-1/2 



(5) 



where k is the eigenvalue of a K 1 \k\ = 
\fj(j + 1) + 1/4 = j + 1/2. Since k can have two 
signs the energy spectrum is two-fold degenerate. In 
equation Q for JL operator the Dirac operator is 
denoted by j. In our approach we follow the simple 
logic: 

- If there is some symmetry that relates two signs of re, 
the corresponding symmetry operator must anticommute 
with K . 

- At the same time, naturally, this operator has to 
commute with the Dirac Hamiltonian. 

Thus the first step will be constructing an operator (s), 
that anticommute with K. In order to find such an oper- 
ator, let generalize the theorem, which was known earlier 
for the Pauli equation. We reformulate it in the following 
form. 

Theorem: Let V be a vector with respect to the angular 
momentum I, i.e., 

[k, Vj] = iSijkVk 
In the vector product form it can be written as 

txV + V x T= 2iV 
Suppose also that this vector is perpendicular to I 



Then K anticommutes with operator \ Y, ■ V j , which is 
scalar with respect to the total J momentum. 



Proof: Let us consider a product ^E -tj ( S ■ V^j . Ex- 
ploiting the known properties of Dirac matrices and con- 
ditions of the theorem, one can establish that 

(t ■ t\ (e • f) = (T> V ) + i (£, fx V ) = 

= i (E, 2iV - V X tj = -2 (f, ■ Vj - i (E, fx V^j . 
Therefore 

(S • f+ ■ v) = - |E • V + i (£,, fx V ) } . (6) 

Now consider the same product in reversed order 
(E • V) (t ■ t\ = (v ■ t\ + i (£, V x I) = i (E, V xtj . 
Hence 

(S • V\(t ■ f+ l) = (E • Vj + i (£, Vxf\ = 
= - (£ • 1 + lVs • V 



In the last step we made use of the equation ©. There- 
fore we have obtained 



(S-Z + l,E-y) =0 (7) 
Now, according to the definition J5J), it follows that 

K (E • V \ = - (£ • V\ K (8) 

Thus the theorem is proved. 

It is evident that the class of anticommuting with K 
(so called, K -odd) operators is not restricted by these 

operators only - any operator of kind O ^E • Vj , where 

O is commuting with K , but otherwise arbitrary, also is 
K -odd. 

It is noteworthy, that the following useful relation 
holds in the framework of constraints of above theorem: 

1 



K 



(Z.v) 



-i(3 E 



J ' 2 



V xl-lxV 



(9) 



which follows from vector's transformation rule men- 
tioned in theorem in the form of vector product. 

One can see that the antisymmetrized vector product 
characteristic of LRL vector appears on the right-hand 
side of this relation. Important special cases, resulting 
from the above theorem include V = f (unit radial vec- 
tor), V = p (linear momentum) and V — A (LRL vector). 
This last one has the following form£ 



A = f - 



2ma 



p x I — I x p 



(10) 



According to there is a relation between above three 
odd operators 

A = Y,-"r+ — [3k(t,-p) (11) 
ma V / 

At this point we are tooled up to derive the hidden sym- 
metry operator, which is done in the following section. 
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IV. THE HIDDEN SYMMETRY OPERATOR 

The second step of our two stage derivation strategy is 
to find K-odd operators that commutes with the Dirac 
Hamiltonian. There remains still considerable freedom 
of choice here because of above mentioned remark about 
operators like O (fl ■ V^j . One can take O into consider- 
ation or ignore it. 

Let us choose 



E • r and K 



(12) 



Because of specific place of LRL vector, it is reasonable 
to assume that E • A - like term will present here. So 
this choice is dictated by ljll|l . Both operators in l|12|l 
are diagonal matrices. After commuting them with non- 
diagonal H we will eventually end up with non-diagonal 
terms. For example, 



Y, ■ r,H 



= -ftKj 5 
r 



(See, exercise 2 below). The resulting matrix on the 
right-hand side is anti-diagonal. Probing the follow- 
ing combination of diagonal and anti-diagonal matrices 
yields: 

A' = x 1 (t-f) + ix 2 K (s ■ + ix 3 Kf 5 f(r) (13) 



Here the coefficients are chosen in such a way, that A' 
be Hermitian, with 27,2,3 are arbitrary real numbers and 
f(r) is an arbitrary scalar function (to be determined 
later). 

The commutator of A' with H is given by 

[A',H] = xi-ftKj 5 - x 2 \k (t-f) - 
r r \ J 

-x 3 f(r)K (s ■ f) - ix 3 2m(3Kj 5 f(r) (14) 



and using this result in the first equation, we obtain that 

1 



x 2 = Xl 

ma 



%3f(r) = — xi 
mr 



(18) 



(19) 



Substitution of flTJl. (fT%|) and (O into U3J) yields the 
following operator that commutes with the Dirac Hamil- 
tonian, 10: 



A' = xi I (s ■ f] - —K (s ■ 0) + —K-v 
\ J ma \ J mr 



(20) 



This operator is a A" -odd and satisfies all conditions of 
the above theorem. 

Thus we constructed the operator, which is associ- 
ated with a hidden symmetry of the Dirac equation in 
Coulomb field. But this operator is not new. One can 
make sure of that (|20|) is a different form of JL operator. 
Indeed, if we turn to usual a matrices making use of def- 
inition Q and taking into account the expression @ for 
the Dirac Hamiltonian, A' can be reduced to the form 
(xi, as unessential common factor, may be dropped) 



A' = 7 5 <^ a ■ r A7 5 (H - (3m) 

' ma 



This expression is nothing but the Johnson-Lippmann 
operator, Q. 



V. PHYSICAL MEANING AND SOME 
APPLICATIONS OF JL OPERATOR 

In order to understand a physical meaning of A' oper- 
ator note that, using JUJl, equation l|20|l may be rewritten 
in the following form: 



Grouping diagonal and antidiagonal matrices separately 
and equating this expression to zero, we obtain equation 



K E-f 



(^2+x 3 f'(r)) + 



+2i(3K^ -xi - mf{r)x 3 = 



(15) 



This equation is satisfied, if diagonal and antidiagonal 
terms become zero separately, i.e., 



A' = t 



2ma 



ft 



p x I — I x p 



AV (21) 

mr 



In the non-relativistic limit, when (3 — > 1, and 7 — > 0, 
this operator reduces to the projection of LRL vector on 
the electron spin direction, A' — » E • A, or because of 
(l ■ Aj = 0, it is a projection on the total momentum, J. 

As for further application, let us calculate the square 
of JL operator. The result is as follows^: 



1 



x i = mf(r)x 3 , 



7X2 



f(r)x 3 (16) 



Integration of the second equation over the interval (r, 00) 
yelds 



X3f(r) = — x 2 
r 



(17) 



A 2 



H 2 



- 1 



(22) 



Because all operators in l|22() commute with each other, 
one can replace them by their eigenvalues. Therefore one 
obtains energy spectrum pure algebraically after specify- 
ing spectrum of A 2 . Since A 2 is positively defined the 
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minimal eigenvalue of A 2 is zero. For this eigenvalue so- 
lution of gives precisely the ground state energy of 
hydrogen atom, 



E = m 1 - 



(Zaf 



1/2 



Full spectrum can be easily derived by well-known ladder 
procedure^ . It leads to the usual Sommerfeld formula 
©■ 

It is worthwhile to note a full analogy with classical me- 
chanics, where closed orbits were derived by calculating 
the square of LRL vector without solving the differential 
equation of motion^. Therefore we are convinced that 
the degeneracy of hydrogen atom spectrum with respect 
to interchange k — > — K is related to the existence of JL 
operator, which in its turn takes its physical origin from 
the Laplace-Runge-Lenz vector. 

It is also remarkable that the same symmetry is re- 
sponsible for absence of the Lamb shift in this problem. 
Inclusion of Lamb shift termsii^ 



A Via 



Aa 2 
3m 2 



In- 



6 3 (r) 



27rm 2 r 3 



(23) 



found by calculating radiative corrections to the pho- 
ton propagator and photon-electron vertex function, into 
the Dirac Hamiltonian breaks commutativity of A with 
H. However it is evident that without radiative correc- 
tions, terms like (|23(l do not appear in the Dirac Hamil- 
tonian, as in the one-electron theory, and as long as only 
Coulomb potential is considered, the appearance of the 
Lamb shift should be always forbidden. 

In conclusion we want to underline once again that 
the hidden symmetry, associated to the Coulomb po- 
tential, governs a wide range of phenomena from plan- 
etary motion till fine and hyperfine structure of atomic 
spectra. Therefore teaching of special aspects of this ex- 
tremely unusual symmetry can extend the capabilities of 
advanced students. 

Finally we present some useful exercises that can help 
students to master calculation of commutators, consid- 
ered in the main text. 



VI. EXERCISES 

Making use of these exercises students will enhance 
their knowledge in technical problems of manipulations 
with Dirac matrices and calculations of some commuta- 
tors. 

Exercise 1. 

Show, that the Dirac's operator commutes with Hamil- 
tonian, i.e. [K, H] — 0, for arbitrary central potential, 
V(r). 



Solution: 

[K,H]= + 1, a-p + /3m + V(r) 

p(p-t\, a ■ p + [f3, a ■ p\ = (3 ("s • l\ (a ■ p) - 

- {a ■ p) (e • t\ + 2f3 (a ■ p) = j3 (e ■ t\ (a ■ p) + 

+(3{d-p) (E-z) + 2j3(a-p) 

perform transitions from alpha to sigma matrices 

[K,H] = /? 7 5 {(e •/) (E-p) + (s -pj (e • tj + 

+2 (s-p)} =/? 7 5 {-2(s-p) -i(z,pxl) + 

+i (f,,pxlj + 2(£-p\\ = 

At the last step we made use of theorem given in the 
text. 

In remaining exercises below calculation of commuta- 
tors are require, that are needed for calculation of com- 
mutator of the symmetry operator with the Dirac Hamil- 
tonian. 

Exercise 2 



Calculate the commutator 
Solution: 

E • r, a • p + /3m + V(r) 



E • r,H 









E • r, a • p 



-,(XjPj 



Let calculate individual terms: 



[Ej, ajPj] — 



-,a jPj 



Ej- [riUjPj] + E; 



, a jPj 



ri = -E ■ a — «E ■ r — ^~ 



3* 5 i 5 

—7 - -7 

r r 



[Ei,ajPj] ^ = 2e ljk a kPj r i = 
= 2ie ljk a k {r iPj - iS tj ) C^J = 
= 2ia k l k f ij = ^ 7 5 (E • t) 

Grouping these terms together, we derive 

2i 



2 1 



7 5 + - 7 5 (E-f) 



= - 7 5 fs.f+ l) =- 1 5 (3K 
r \ I r 



Exercise 3 

Calculate the commutator 



K 



(e-p), 



H 
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Solution: 

K(Z-p),H 



K 



[K,H] (E-p) 



K 



E • p,a ■ p + (3m, + V(r) 



K 



e • p,d ■ p 



-mK 



K 



E • p, V{r) 



Let calculate individual commutators 



E • p, a ■ p 



[Ei, on] PiPj = 2i ", ,/ '>/,/»,//, = 



Exercise 4 Calculate commutator [i^7 5 /(r),H] 
Solution: 



[X 7 5 /M,^] - K [j 5 f(r),H] + [K,H} 7 5 f(r) = 
= K~f> [/(r), H] + K [ 7 5 , H] f(r) = [f(r), a ■ p\ + 
+K ■ 2m 7 5 ^/(r) = K^ia. ■ ~rf(r) + 2mif 7 5 /3f(r) = 
= iK fE • /'(r) + 2mK~/ 5 Pf(r) 



E-p,/3 =0 



E • £ y(r) = E, : [p t , V(r)} = -t^V'(r) 



■i (E • r \ V'(r) 



For the Coulomb potential V(r) = —a/r — > V'(r) 
a/r 3 and hence 



K 



(2-*H=-£(S.*) 



These relations are applied above in calculating of 
equation l|14() . 
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